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ABSTRACT 



In this paper, we investigate the effect of cooling on the X-ray properties of galaxy 
clusters. We have performed N-body, hydrodynamical simulations both with and with- 
out the effects of radiative cooling, but neglecting the effects of star formation and 
feedback. We show that radiative cooling produces an inflow of high-entropy gas from 
the outer parts of the cluster, thus raising the cluster temperature and decreasing the 
X-ray luminosity. With radiative cooling clusters are on average three to five times 
less luminous in X-rays than the same cluster simulated without cooling. However, 
we do not produce a large constant-density core in either the gas or the dark matter 
distributions. Our results contradict previous work in which cooling raises the X-ray 
luminosity and deposits an unreasonably large amount of mass in the central cluster 
galaxy. We achieve this by selecting our numerical resolution in such a way that a 
reasonable fraction of the baryonic material cools and by decoupling the hot and cold 
gas in our simulations, a first step towards modelling multiphase gas. We emphasise 
that globally cooling a sensible amount of material is vital and the presence or absence 
of massive central concentrations of cold baryonic material has a dramatic effect upon 
the resultant X-ray properties of the clusters. 
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1 INTRODUCTION 

Clusters of galaxies are the largest virialised structures in 
the Universe, evolving rapidly at recent times in many pop- 
ular cosmological models. Even at moderate redshifts the 
number of large dark matter halos in a cold dark matter 
Universe with a significant, positive cosmological constant 
is higher than in a standard cold dark matter Universe and 
it is precisely because both the number density and size of 
large dark matter halos evolve at different rates in popu- 
lar cosmological models that observations of galaxy clusters 
provide an important discriminator between rival cosmolo- 
gies. 

The advent of X-ray satellites opened up a whole new 
area in observational astronomy. Hot gas, typically at tem- 
peratures of 10 7-8 K sitting in the deep potential wells of 
galaxy clusters emits radiation via thermal bremsstrahlung. 
This emission is heavily biased towards the central regions 
of the cluster because the flux is weighted as the gas density 
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squared. Given the gas temperature and the X-ray surface 
brightness, the gas column density and a spherically sym- 
metric gas density profile can be estimated. If the hot gas 
is assumed to reside close to hydrostatic equilibrium within 
the dark matter halo the underlying dark matter density 
distribution can be derived. 

Since the early work of Abramopoulos & Ku (1983) and 
Jones & Forman (1984) on the radial density profiles of large 
galaxy clusters, debate has raged about the presence or ab- 
sence of large, ~ 250 kpc, constant density cores in the 
X-ray emitting gas. Unfortunately X-ray imaging is notori- 
ously difficult because of the inherent large beam size. In 
addition, the centre of the emission region is not easy to 
determine and any centering error also acts to smooth out 
any central increase when the results are azimuthally aver- 
aged and plotted as a radial profile (Beers & Tonry 1986). 
More recent high resolution images of several clusters have 
helped to resolve this controversy. Some galaxy clusters do 
indeed appear to exhibit a large, resolved core but others 
have a much smaller core, close to the resolution threshold of 
the instrument. Both White, Jones & Forman (1997) (who 
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studied 207 clusters imaged by the Einstein observatory) 
and Peres et al. (1998) (who looked at ROSAT observations 
of the flux-limited sample of clusters provided by Edge et 
al. 1990) suggest that this dichotomy relies on the presence 
or absence of a cooling flow. Clusters with a cooling flow ap- 
pear to have small cores (around 50/i _1 kpc) whilst clusters 
without cooling flows have much larger cores. The work of 
Allen (1998) on the discrepancy between the large X-ray core 
radii and the small core radii deduced from strong lensing 
observations (e.g. Kneib et al. 1996) also reached the con- 
clusion that cooling flow type clusters had small core radii 
in their matter distributions. Because the X-ray flux rises as 
the local density squared, the total X-ray emission from a 
cluster is very sensitive to the central density. 

High resolution N-body simulations of galaxy clusters 
(Moore et al. 1998) produce a radial dark matter profile that 
has no core. The radial profile continues to rise until the res- 
olution threshold is reached, well within the required radius 
if the gas is to trace the dark matter and still reproduce the 
X-ray observations. The production of a central constant 
density dark matter core has been a long standing problem 
for collisionless dark matter simulations (Pearce, Thomas & 
Couchman 1993, Navarro & White 1994), although in pre- 
vious work the resolution threshold was still close to the 
observed core sizes and so until the latest high resolution 
studies this was still a tentative result. 

Recently several groups have tried to reconstruct the 
radial temperature profile of galaxy clusters. Markevitch et 
al. (1998) used ASCA data from 30 clusters and concluded 
that the temperature falls steeply at large radii. However, 
Irwin, Bregman & Evrard (1999) analysed ROSAT PSPC 
images of 26 clusters and concluded that the radial temper- 
ature profiles were generally flat out to the virial radius. 

The inclusion of a gaseous component into simulations 
allows the previously assumed relationship between the gas 
and the dark matter to be derived directly. The first study 
of this type (Evrard 1990) was carried out without the ef- 
fects of radiative cooling but reproduced well many theo- 
retical predictions such as a bias between the dark matter 
and the baryonic material. Similar non-cooling simulations 
have proved popular (Cen & Ostriker 1994, Kang et al. 1994, 
Bryan et al. 1994, Navarro, Frenk & White 1995, Bartel- 
mann & Steinmetz 1996, Eke, Navarro & Frenk 1998, Bryan 
& Norman 1998) and such a simulation formed the basis 
of the Santa Barbara project in which a dozen groups sim- 
ulated the formation of the same galaxy cluster (Frenk et 
al. 1999). 

Extending this work to include a dissipative component 
has proved difficult because the formation of galaxies intro- 
duces many additional physical effects. Metzler & Evrard 
(1994) and Evrard, Metzler & Navarro (1996) circumvented 
this by introducing a galactic component by hand into a 
simulation that didn't follow radiative cooling of the gas to 
study the effects of feedback on the cluster profile. Fully 
consistent attempts to follow the radiative cooling of the 
hot intra-cluster gas have only recently been achievable be- 
cause of the extra computational overhead involved. Early 
attempts to include the effects of radiative cooling (Thomas 
& Couchman 1992, Katz & White 1993, Evrard, Summers 
& Davis 1994, Frenk et al. 1996) either suffered from poor 
resolution, focussed on the galactic population or suffered 
from overmerging effects. 



In a now classic paper Katz & White (1993) examined 
the effect of radiative cooling on the X-ray profile of a sin- 
gle galaxy cluster, a study repeated recently by Lewis et 
al. (1999). Their simulated cluster has properties that are 
not observed; a Virgo sized cluster with a super-massive cen- 
tral galaxy and an enormous 400 M Q / yr associated cooling 
flow. As they point out, massive brightest cluster galaxies 
of this size are observed within the Universe, as are massive 
cooling flows, however they fail to stress their rarity, partic- 
ularly in objects of similar size to the Virgo cluster. More 
recently Suginohara & Ostriker (1998) simulated a different 
cluster and also produced an object which has properties 
they themselves admit are unobserved - "The high resolu- 
tion simulation resulted in a gas density profile steeply rising 
toward the center, with consequent very high X-ray luminos- 
ity; however, these properties are not observed". They sug- 
gest that feedback of energy from supernovae might account 
for the discrepancy. In this paper, we obtain quite different 
results because, unlike previous studies, a realistic fraction 
of the baryonic material has cooled to form galaxies. Large, 
bright central cluster galaxies can have a dramatic effect on 
the cluster potential and consequently the X-ray properties. 
A reasonable treatment of these objects is extremely impor- 
tant in studies of this type. To prevent too much material 
cooling to form the central object within each halo we em- 
ploy a modified form of smoothed-particle hydrodynamics 
— this will be described further in Section 2, below. 

If the effects of cooling are included in the models then 
one might expect the entropy of the intracluster medium to 
decrease. Paradoxically, several mechanisms have been sug- 
gested that may produce large, constant-density gas cores by 
raising the entropy of the gas at the centre of the cluster: 

• Radiative cooling is very efficient in small dark matter 
halos because the cooling time is less than the dynamical 
time (White & Rees 1978). These knots of cold, dense gas 
can be equated with proto-galaxies and, as they are dense 
and collapsing, they may be reasonably expected to produce 
stars (Katz 1992). Star formation leads to energy feedback 
into the interstellar medium via supernovae explosions (Katz 
1992, Mihos & Hernquist 1994, Navarro & White 1994, Ger- 
ritsen & Icke 1997). Unless the gas immediately recools, this 
heating acts to increase the entropy of the surrounding ma- 
terial, pushing it onto a higher adiabat and preventing it set- 
tling to the very high densities and temperatures required 
for it to trace the underlying dark matter (Wu, Fabian & 
Nulsen 1998). 

• The presence of galaxies orbiting within the cluster po- 
tential acts to stir up the gas, heating it as friction and 
turbulence dissipate the galaxies' velocity, simultaneously 
producing velocity and spatial bias in the galaxy distribu- 
tion (Frenk et al. 1996). This effect is most pronounced in 
the centre of the cluster. 

• A third mechanism for producing a core is that radia- 
tive cooling of the gas at the centre of each potential well 
acts as a drain on the low entropy material (which cools pref- 
erentially). If the remaining gas cannot cool rapidly enough, 
a core would develop because only high entropy material 
remains, an effect postulated in the first paper to include 
radiative cooling (Thomas & Couchman 1992) and later re- 
iterated by Waxman & Miralda-Escude (1995) and Bower 
(1997). It is this mechanism that we investigate in this pa- 
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per. We show that the entropy of the intracluster medium is 
indeed increased, that this leads to a greatly reduced X-ray 
luminosity, but that it does not give a large, constant-density 
core. 

The remainder of this paper is laid out as follows: in 
Section 2 we present the large hydrodynamical simulations, 
both with and without the effects of radiative cooling, that 
we have performed; in Section 3 we extract radial density 
and temperature profiles for the 20 largest galaxy clusters 
within each simulation and contrast the profiles with the 
underlying dark matter distribution; this is followed in Sec- 
tion 4 by a discussion of our findings. 



2 THE SIMULATIONS 

The simulations that we have carried out use the adaptive 
particle-particle, particle-mesh (AP 3 M) method (Couchman 
1991) coupled to the smoothed particle hydrodynamics 
(SPH) technique (Gingold & Monaghan 1977, Lucy 1977) 
to follow 2 million gas and 2 million dark matter particles in 
a box of side 100 Mpc (Couchman, Thomas & Pearce 1995, 
Pearce & Couchman 1997). 

We have performed simulations in two types of flat cold 
dark matter cosmology, one standard (SCDM) and one with 
a cosmological constant ( ACDM) , with the same parameters 
assumed by Jenkins et al. (1998) (Q = 1.0, A = 0.0, h=0.5, 
a s = 0.6 for the former and SI = 0.3, A = 0.7, h=0.7, 
o~s = 0.9 for the latter). The baryon fraction was set from 
Big Bang nucleosynthesis constraints, Slbh 2 = 0.015 (Copi, 
Schramm & Turner 1995) and we have assumed an unevolv- 
ing gas metallicity of 0.3 times the solar value. These pa- 
rameters produce a gas mass per particle of 2 x 10 9 Mq in 
each case and are summarised in Table 1. The dark mat- 
ter mass is only slightly lower than that given by Steinmetz 
& White (1997 — their Equation 9, but note that there is a 
typographical error) at which artificial 2-body heating bal- 
ances radiative cooling. Thus we expect there to be some 
numerical heating in our simulations. However, as we have 
chosen to neglect real heat sources, such as supernovae in 
galaxies, this is of little importance and does not affect our 
conclusions which concern the differences between runs with 
and without radiative cooling. 

Since we smooth over 32 SPH particles, the smallest 
gaseous object that can be effectively resolved has a mass of 
6.4 x 10 10 M . We employ a comoving /3-spline gravitational 
softening equivalent to a Plummer softening of 10/i -1 kpc 
for redshifts < z < 1.5 (2.5 for ACDM) — at earlier times 
the softening has a fixed physical size, with the minimum 
SPH resolution set to match this. We note that a spatial 
resolution of 10/i -1 kpc is over twice the typical scale length 
of elliptical galaxies and this may lead to enhanced tidal 
disruption, drag and merging within the largest clusters of 
objects. However, the force softening cannot be reduced fur- 
ther without introducing 2-body effects. A smaller soften- 
ing would also lead to a further increase in the number of 
timesteps required; we already require around 10000 for each 
cooling run. 

In addition to these two simulations which both in- 
cluded the effects of radiative cooling, we repeated the 
ACDM model without cooling. This simulation has input 
parameters close to those used by Eke, Navarro & Frenk 
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Table 1. The main parameters for each cosmology. The param- 
eters for the cooling and non-cooling runs were identical. 

(1998) who used the same cosmology and the resultant clus- 
ters look very similar and exhibit similar X-ray properties 
(see Fig. for a comparison). 

The properties of the galaxies in the two simula- 
tions with radiative cooling have been described in Pearce 
et al. (1999, 2000). They clearly demonstrate an accept- 
able match to both the spatial and luminosity distribu- 
tion of observed galaxies. This was achieved by employ- 
ing three numerical approximations: a mass-resolution of 
6.4 x 10 10 Mq below which objects cannot cool efficiently, a 
length-resolution (ie softening) of Wh kpc, and decoupling 
of the hot, halo gas from the cold galactic gas. Improving 
the mass and/or length resolution would increase the frac- 
tion of cold gas in our simulations, producing galaxies that 
were too luminous to match the observations. This would 
then necessitate the introduction of feedback mechanisms 
that would over-complicate our model. 

Decoupling the hot halo gas is a new innovation that 
we feel vastly improves the ability of SPH to model fluids 
in which there are large density contrasts. Without it, hot 
gas particles have their density overestimated in the vicinity 
of cold gas and too much material cools to form the central 
galaxy. Normal cooling of the intracluster medium at tem- 
peratures above 10 K is still handled correctly and galaxy- 
galaxy mergers and viscous drag on the galaxies as they 
orbit within the halo are retained. As Pearce et al. (1999) 
show, such a procedure produces a set of galaxies that fit 
the local K-band number counts of Gardner et al. (1997). 
The brightest cluster galaxies contained within the largest 
halos are not excessively luminous for a volume of this size, 
unlike those found in previous work (Katz & White 1993, 
Lewis et al. 1999). The fraction of the baryonic material that 
cools into galaxies within the virial radius of the large halos 
in our simulation is listed in Tables 3 & 4 and is typically 
around 20 percent. This is much less than the unphysically 
high value of 40 percent reported by Katz & White (1993). 
Decoupling of the hot phase produces a galactic population 
and cold gas fraction that is well matched to the observa- 
tions. 



3 RESULTS 

3.1 Extracting objects 

For the purposes of this paper we are interested in only the 
largest objects within each simulation, as only these contain 
sufficient mass to produce the deep potential wells required 
to retain hot, X-ray emitting gas. We centre each cluster on 
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the peak of the hot gas density, a position that coincides 
with the centre of the X-ray emission. This prevents the 
introduction of an artificial constant-density core which may 
arise with any other choice of centre — for clusters with 
significant substructure the centre-of-mass can lie a long way 
from the centre of the X-ray emission. 

The virial radius for each of our clusters was defined as 
the spherical region, surrounding each cluster centre, that 
enclosed an overdensity of 178 for SCDM and 324 for ACDM 
(Eke, Cole & Frenk 1996). Each catalogue was then cleaned 
by ordering it in size and deleting the smaller of any over- 
lapping clusters. The 20 most massive clusters in each of the 
catalogues were then used for the work presented here. The 
properties of the clusters are presented in Tables 2-4. The 
16 largest clusters found in the non-cooling simulation are 
found in the list of the 20 largest clusters in the other two 
models. The index of the matching cluster from the non- 
cooling run is given in Tables 3 & 4. 

Each of the extracted clusters was checked for substruc- 
ture by comparing the centre of the X-ray emission to the 
median position of the particles within the virial radius, a 
statistic that has been shown to be a useful indicator of the 
presence of substructure by Thomas et al. (1998). The re- 
sults of this test are shown in Tables 2-4. All the clusters 
with an offset of more than 7 percent of the virial radius (6 
in each case) were noted and are shown as dotted lines on 
Figures | & |^ and as open symbols on Figure (jj. 

3.2 Dark matter density profiles 

The radial dark matter density profiles for the 20 largest 
objects in each cosmology are shown in Figure [jj Shown 
is the mean density within spherical shells, the innermost 
shell plotted in the Figure containing at least 64 particles. 
The dark matter profiles of those clusters without significant 
substructure are similar within each cosmology. In the non- 
cooling run, the innermost bin of each density profile shows 
a flattening. This is a resolution effect: the radial extent of 
the bin has to be very large in order to accommodate 64 
particles. 

Both gas and dark matter density profiles for the largest 
object in each of the runs are shown in Figure ^. The dark 
matter profile for the ACDM run without cooling is reason- 
ably well fit by an NFW profile (Navarro, Frenk & White 
1997). With cooling the dark matter density profile of the 
clusters is not well fit by the NFW formula, because the 
asymptotic slope in the central regions is steeper than -1. 
This is because, once cooling has been implemented, the 
large galaxy that forms at the centre of each cluster acts 
to draw in more dark matter and steepen the profile signifi- 
cantly in the inner regions. A similar effect is seen for most of 
the other clusters, although in the cooling run several show 
a drop in density in the innermost bin. This indicates that 
the peak X-ray emissivity sometimes comes from a galaxy 
that is not located at the centre of the cluster. 



3.3 Gas entropy profiles 

Gas entropy profiles (and also density and temperature pro- 
files, below) were obtained using only those particles with a 
temperature exceeding 12 000 K. Typical gas temperatures 



exceed 10 7 K for halos in this mass range and we wish to ex- 
clude cold gas which lies within galaxies or recently tidally 
disrupted objects. If the cold material were included, there 
would be a large density spike at the centre of each of the 
clusters (which all have a central galaxy). This object does 
not contribute to the X-ray emission because the material 
it contains is very cold compared to the surrounding hot 
halo (although the increased depth of the local potential 
can help to confine dense, hot gas which can affect the bolo- 
metric X-ray emission) . The specific entropy profile is shown 
in Figure |. We plot the quantity (T/ K)/(p/p) 2/3 , where T 
is the temperature and p the density, measured in units of 
the mean gas density, p. 

Let us contrast the results for the ACDM runs with and 
without cooling. Firstly, note that the entropy at the virial 
radii is very similar in each case — this is because cooling has 
had little effect at these large radii. Between the virial radius 
and about 0.2 times the virial radius (less for the largest 
cluster) , the entropy profiles for the cooling run are shallower 
than in the non-cooling run. This confirms the hypothesis 
of Thomas & Couchman (1992) that cooling is able to raise 
the entropy of the intracluster medium by dragging in high- 
entropy material from the outer regions of the cluster. 

Within about 0.2 virial radii, the entropy profiles again 
steepen — it is within this "cooling radius" that the cooling 
time is short enough to allow significant cooling of the gas 
within the lifetime of the cluster. By the time we get to the 
innermost bins in the Figure, there seems to be a spread 
in the entropy of the clusters in the cooling run, with some 
having higher entropy and some lower entropy than the cor- 
responding clusters in the non-cooling run. 

The SCDM run with cooling exhibits similar entropy 
profiles to the ACDM run with cooling. 

3.4 Gas density profiles 

The radial gas density profiles are displayed in Figure ^. 
Without cooling these profiles look very similar to those ob- 
tained by Eke et al. (1998). Within about 0.1 virial radii, the 
gas profiles are significantly shallower than the correspond- 
ing dark matter profiles of Figure |l]. This is indicative of 
the fact that the hot gas has a higher specific energy than 
the dark matter and, as previous authors have found {e.g. 
Navarro et al. 1995, Eke et al. 1998), there is more dark mat- 
ter than gas (relative to the cosmic mean) within the virial 
radius. There is a general tendency for all the profiles to 
flatten considerably in the innermost bin. Once again, this 
is due to inadequate resolution and we do not attach any 
significance to it. 

The effect of cooling is to lower the gas density at all 
radii within the virial radius. The suppression is greatest, a 
factor of three, at about 0.1 times the virial radius, roughly 
corresponding to the kink in the entropy profiles seen in 
Figure ^. Although the density gradients are shallower, they 
do not roll over into constant-density inner core regions. In 
fact, for the larger clusters, the density continues to rise 
further into the centre of the cluster than before, so that 
the central density is close to that in the non-cooling case. 

The profile of the largest object in both gas and dark 
matter for each of the runs is shown in Figure PJ. The arrow 
indicates a radius of 100/i _1 kpc. Without cooling, the gas 
density is shallower than that of the dark matter within 0.1 
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Table 2. Data for the 20 largest clusters (by mass) in the ACDM run without cooling. Quoted are the cluster number, the virial mass in 
units of 10 14 /i _1 Mq, the virial radius in Mpc, the value of the substructure parameter defined in the text (a larger number implies 
more substructure), the number of dark matter particles within the virial radius, the number of gas particles within this radius, the 
emission weighted mean temperature (in keV/fc) and the bolometric X-ray luminosity in units of 10 43 erg/s/h 2 
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Table 3. Data for the 20 largest clusters (by mass) in the ACDM run which includes cooling. Quoted are the cluster number, the number 
of the matching cluster in the ACDM run without cooling (if one exists), the virial mass in units of 10 14 /i — 'Mq, the virial radius in 
h~ 1 Mpc, the value of the substructure parameter defined in the text, the number of dark matter particles within the virial radius, the 
number of gas particles above and below 12000 K within this radius, the number of cold particles within 100/i~ x kpc of the centre, the 
fraction of the baryonic material within the virial radius that has been able to cool, the emission weighted mean temperature (in keV/fc) 
and the bolometric X-ray luminosity in units of 10 43 ergs/s//i 2 . 

times the virial radius, but this inner, resolved slope of the 
density profile is still p oc r" 1 with no sign of a constant- 
density core. As the temperature is approximately constant 
within this region (see Figure |B|), the X-ray luminosity is con- 
vergent and dominated by emission from around 200 h~ kpc 
(0.1 times the virial radius). 

With cooling, the largest cluster exhibits a central den- 
sity spike due to the presence of a massive central galaxy. 
This hot gas has a very steep radial density profile, p oc r -3 , 



and would be classified observationally as a cooling flow of 
60h~ 2 Mq/yi onto the central cluster galaxy. Between radii 
of about 40/i _1 kpc and l/i _1 Mpc, the density profile is a 
power law, p tx r" 1 ' 4 , steepening at larger radii. Thus the 
X-ray luminosity (excluding the cooling flow) comes from a 
much more extended region than in the non-cooling case. 

In conclusion, the gas density has been reduced by the 
influx of high-entropy material, as expected. However, this 
has not given rise to constant-density inner cores. In fact, if 
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Table 4. Data for the 20 largest clusters (by mass) in the SCDM run which includes cooling. Quoted are the cluster number, the number 
of the matching cluster in the ACDM run without cooling (if one exists), the virial mass in units of 10 14 /i — 1 Mq, the virial radius in 
h~ 1 Mpc, the value of the substructure parameter defined in the text, the number of dark matter particles within the virial radius, the 
number of gas particles above and below 12000 K within this radius, the number of cold particles within lOOh^ 1 kpc of the centre, the 
fraction of the baryonic material within the virial radius that has been able to cool, the emission weighted mean temperature (in keV/fc) 
and the bolometric X-ray luminosity in units of 10 43 ergs/s//i 2 . 



anything, the density profiles now continue as a power-law 
closer into the centre of the clusters. 

3.5 Radial temperature profiles 

Radial temperature profiles are shown in Figure [|. The tem- 
perature profiles for the relaxed clusters from the ACDM 
run without cooling are typical of those found in previous 
work (see for example Eke et al. 1998 and references therein). 
They rise inwards from the virial radius by about a factor of 
two, peaking at about 0.1 times the virial radius and then 
declining again, very gradually, in the cluster centre. 

Cooling makes little difference to the temperature pro- 
files, except that corresponding clusters in the ACDM runs 
reach a higher peak temperature when cooling is imple- 
mented, due to the inflow of higher entropy gas. The temper- 
atures are very similar at the virial radius, but are about 1.5 
times higher at their peak than before. Two clusters show 
a precipitous decline in temperature in the cluster centre, 
one of these being the largest cluster — this is evidence for a 
cooling flow. 

The SCDM results are very similar to those for ACDM. 

3.6 X-ray luminosity profiles 

We follow Navarro et al. (1995) in using the following esti- 
mator for the bolometric X-ray luminosity of a cluster, 

**=4X10»X; (f)(|)' ergs- (1) 

where the density is in units of overdensity (relative to the 
mean gas density in the box, 2.86 x 10 -31 g/cm 3 ) and the sum 
extends over all the gas particles with temperatures above 
12000 K. 



The total X-ray luminosity within the virial radius of 
each of the clusters is listed in Tables 2-4. For the simulation 
with cooling the clusters are several times less luminous than 
those from the corresponding non-cooling run. This contra- 
dicts the previous results of Katz & White (1993), Sugi- 
nohara & Ostriker (1998) and Lewis et al. (1999) who all 
found the X-ray luminosity increased if cooling was turned 
on. The reason for the discrepancy is, once again, the fact 
that we have decoupled the hot and cold gas, thus greatly 
suppressing the cooling of the inflowing, high-entropy gas in 
our simulations compared to previous ones. This causes a 
large reduction in the mass of the brightest cluster galaxy 
compared to those produced by previous work. Our galaxies 
have reasonable luminosities, mass-to-light ratios and num- 
ber counts for a volume of this size. 

Note that estimates of X-ray luminosity from the non- 
cooling run are not really meaningful. A radiation rate of 
this magnitude can only be sustained for a short time before 
depleting the intracluster medium of gas. The cooling runs 
produce a more physically self-consistent X-ray luminosity 
because the radiative effects are taken into account. The X- 
ray luminosity within the cooling radius is approximately 
equal to the enthalpy of the gas divided by the age of the 
cluster. 

We plot these bolometric luminosities as a function of 
radius for each of our clusters in Figure ^| Clearly the rela- 
tive contribution to the total X-ray emission from different 
radii is very different for the cooling and non-cooling simu- 
lations. Without cooling all the relaxed clusters show very 
similar emission profiles, with only a small contribution to 
the total emission coming from the very centre. These pro- 
files are mostly well resolved, as claimed by Eke et al. (1998) 
for simulations of clusters with this particle number. 

Once radiative cooling is turned on the radial emission 
profiles span a much broader range. For two of the clus- 
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Fraction of virial radius 

Figure 1. The radial dark matter density profiles of the 20 largest, distinct dark matter halos extracted from each of the three simulations. 
Plotted is the mean density within successive spherical shells. Each profile has been scaled to the virial radius and overdensity relative 
to the cosmic mean for the dark matter. Those profiles marked as dotted lines are for clusters that contain significant substructure (see 
text). All the profiles are started from the point where 32 particles are enclosed. 



ters, a central cooling flow type emission is clearly visible — 
contributing 50 percent and 80 percent of the total X-ray 
flux. Although the cooling flow region is not well-resolved, 
it cannot be much larger without depleting the intracluster 
medium of even more gas. For each of the other clusters, 
the radius enclosing half of the total emission is much larger 
than that for the simulation without cooling. 

3.7 Lx-Tx relation 

There has been much debate in the literature centering on 
the X-ray cluster Lx versus Tx correlation. The emission 
weighted mean temperature is plotted against the bolomet- 
ric luminosity within the virial radius for all our clusters 
in Figure |^. The filled symbols represent the relaxed clus- 
ters and the open symbols denote those clusters that show 
significant substructure. 



The effect of cooling is, in general, to slightly raise the 
temperature but to greatly reduce the X-ray luminosity. Ex- 
ceptions are the cooling flow clusters where the large amount 
of emission from gas cooling onto the central galaxy gives 
rise to a lower temperature than in the non-cooling case. 
There are two of these, easily visible on the plot, in the 
ACDM run; we have plotted their new locations, when the 
central cooling flow is omitted, using crosses linked to the 
old location via arrows. In both cases less than 2 percent 
of the hot, X-ray emitting particles were excised to make 
this calculation (606 for cluster Al and 20 for cluster A5). 
In the case of A5 we have caught a transient event - a small 
amount of gas has been reheated by a merger between a 
satellite and the brightest cluster galaxy. This gas is in the 
process of rapidly cooling back onto the central object, emit- 
ting large amounts of X-rays. 

All 3 sets of clusters display a positive correlation be- 
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Figure 2. The radial dark matter (solid line) and gas (dotted line) density profile within the virial radius of the largest dark matter halo 
extracted from each simulation. Each curve is scaled to the virial radius of the halo and overdensity of the appropriate phase relative to 
the mean cosmic density of that phase. In all three runs there is less gas within the virial radius than dark matter, relative to the cosmic 
mean of each species. No constant density core is found in the cooling runs. The arrow indicates a scale of 100/i _1 kpc. 



tween Lx and Tx, although there are insufficient numbers 
to tie the trend down very tightly. It is clear from the com- 
ments in the preceding paragraph that the nature of the 
correlation depends critically upon whether one removes the 
cooling flow emission or not. We believe that a clearer pic- 
ture arises if this is done. 

The regression line in Figure [?] is from Eke et al. (1998) 
and corresponds to Lx oc T x . Our non-cooling clusters fit 
reasonably well with this relation. The cooling clusters lie 
below this line. Given that we expect cooling to be less im- 
portant in the most massive clusters (the absolute value of 
the cooling time and the ratio of the cooling time to the 
dynamical time both increase with cluster mass), then we 
would expect the clusters with radiative cooling to lie closer 
to the regression line at higher Tx ■ Thus the effect of cooling 
should be to steepen the Lx~Tx relation. We hope to test 
this with more simulations of higher mass clusters. 



Also plotted in Figure (?) are the observed data from 
David et al. (1995). Our clusters are smaller and cooler be- 
cause they are not very massive (due to our relatively small 
computational volume), thus it is hard to assess whether the 
cooling or the non-cooling clusters give a better fit to the 
data. SCDM clusters, at least with the simulation parame- 
ters we have chosen in this paper, do not seem to provide 
such a good match to the data. 



4 CONCLUSIONS AND DISCUSSION 

We have performed two N-body plus hydrodynamics simula- 
tions of structure formation within a volume of side 100 Mpc, 
including the effects of radiative cooling but neglecting star 
formation and feedback. By repeating one of the simulations 
without radiative cooling of the gas, we can both compare to 
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Figure 3. The radial specific entropy profiles of the 20 largest, distinct halos extracted from each of the three simulations. Plotted is 
the mean value of T/n 2 ' 3 within successive spherical shells. Each profile has been scaled to the virial radius. Those profiles marked as 
dotted lines are for clusters that contain significant substructure (see text). 



previous work and study the changes caused by the cooling. 
A summary of our conclusions follows. 

(a) Without cooling our clusters closely resemble those 
found by previous authors (Eke et al. 1998 and references 
therein), with dark matter density profiles that closely fol- 
low the universal formula proposed by Navarro et al. (1995). 

(b) With cooling, the formation of a central galaxy within 
each halo acts to steepen the dark matter profile, supporting 
the conclusion of the lensing studies (e.g. Kneib et al. 1996) 
that the underlying potential that forms the lens has a small 
core. This galaxy may not be located exactly at the cen- 
tre of the X-ray emission, sometimes being offset by up to 
SO/i^kpc. The inner slope of the density profile is then 
steeper than that suggested by NFW, closer to that found 
by Moore et al. (1998) from high resolution N-body simula- 
tions (but note that these did not include cooling gas) . 

(c) We confirm the results of Eke et al. (1998) (and previous 
studies quoted therein) that without radiative cooling the 



gas density profile turns over at small radii, that the radial 
temperature profile drops by a factor of two between its 
peak value and that obtained at the virial radius and that 
the baryon fraction within the virial radius is lower than the 
cosmic mean. 

(d) Cooling acts to remove low-entropy gas from near the 
cluster centre, triggering the inflow of higher entropy mate- 
rial. The entropy excess compared with the non-cooling run 
is greatest at about 0.2 times the virial radius, because ra- 
diative cooling lowers the entropy of the gas near the centre 
of the cluster. 

(e) We stress the importance of correctly modelling the 
central cluster galaxy. The resultant X-ray properties of 
the cluster are very dependant upon the centre and if too 
much material cools into the base of the potential well large 
amounts of hot, dense gas can be confined, producing enor- 
mous X-ray fluxes. We have specifically tailored our models 
to both globally cool a reasonable fraction of material and to 
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Figure 4. The radial gas density profile of the 20 largest halos found within each simulation. Those halos containing significant 
substructure are shown as dotted lines. Plotted is the mean gas density within successive spherical shells. All the halos are scaled 
to the virial radius and gas overdensity relative to the baryonic cosmic mean. 



circumvent problems encountered by SPH when faced with 
large density jumps which can lead to high rates of gas cool- 
ing onto the central objects. 

(f) In cooling clusters, the gas density is reduced, by a maxi- 
mum of about a factor of three at 0.1-0.2 times the virial ra- 
dius. The density profile more closely resembles a power-law 
than in the non-cooling run. A few clusters show a central 
density spike, indicative of a cooling flow onto the central 
cluster galaxy (e.g. Fabian, Nulsen & Canizares 1991). 

(g) The temperatures of the cooling clusters show a sig- 
nificant fall between the point where the peak values are 
obtained (around 0.1 virial radii) and the virial radius. The 
observational evidence is somewhat divided here. Our results 
are in agreement with Markevitch et al. (1998) who find that 
the temperature profile of galaxy clusters are steeply falling. 
However, Irwin et al. (1999) recover isothermal temperature 
profiles out to the virial radius from an averaged sample of 
26 ROSAT clusters. It is important to try to clear up this 



observational controversy as isothermal temperature profiles 
are not seen in our models but are often used for X-ray mass 
estimates, measurements of Q b (and hence £2) and theoreti- 
cal arguments for the Lx~Tx relation. 

(h) Cooling acts so as to increase both the mass-weighted 
and observed X-ray temperatures of clusters. The peak tem- 
peratures are raised by a factor of about 1.5 and the tem- 
perature gradient between the peak and the virial radius is 
correspondingly increased. 

(i) The bolometric luminosity for the clusters with radiative 
cooling is around 3-5 times lower than for matching clusters 
without it. Except for the cooling flow clusters, the X-ray 
luminosity profile is less centrally concentrated than in the 
non-cooling case with a greater contribution coming from 
larger radii. 

(j) Cooling flow clusters are easily distinguished from non- 
cooling flow clusters in the Lx~Tx plane. The former are 
more luminous and cooler than the latter (Fabian et al 1994, 
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Figure 5. The radial gas temperature profile found within each of the 20 halos extracted from each cosmology. Clusters with significant 
substructure are marked as dotted lines. All the halos are scaled to the virial radius and temperatures (calculated as the mass-weighted 
mean temperature of the gas within successive spherical shells) are in Kelvin. With cooling the cluster temperature reaches a slightly 
higher maximum at a larger fraction of the virial radius whereas the temperature at the virial radius is very similar. 



Allen & Fabian 1998). Some of this difference results from 
X-ray analysis methods (Allen & Fabian 1998) but may also 
be caused by actual physical differences in the mass distri- 
bution of clusters. We suggest that, while interpretation of 
this relation would be made simpler if the cooling flow were 
excluded before determining the X-ray properties of clusters, 
this should be done with caution. 

(k) The clusters from the non-cooling simulation lie on the 
Lx oc T x regression line of Eke et al. (1988), whereas those 
from the cooling run lie below it. We suggest that, as cooling 
is likely to be less important in more massive clusters, the 
effect of cooling will be to steepen the relation. This remains 
to be tested with simulations of more massive clusters. 
(1) The very large core radii (~ 250h~ 1 kpc) observed in 
some clusters are not seen in our simulations. It is possible 
that such events are rare, or occur preferentially in massive 
clusters, in which case our cluster sample may simply be too 



small to contain such objects. Smaller core radii of around 
~ 50 h- 1 kpc are close to our resolution limit and so are not 
ruled out by our results, although we find no evidence for a 
constant-density region in the centre of any of our clusters. 

At this point, we should remind the reader that our 
simulations are designed merely to investigate the effect of 
radiative cooling on the X-ray properties of the intracluster 
medium and do not set out to explore the complete range 
of physical processes going on in clusters. In particular, we 
ignore the possibility of energy injection into the intraclus- 
ter medium. Early star formation at high redshift and the 
subsequent energy feedback from supernovae could act to 
preheat the gas which falls into the potential wells of galaxy 
clusters, effectively raising its entropy to such a level that it 
cannot reside at the high densities required to trace the dark 
matter into the base of the potential well, and giving rise to 
the constant-density cores that our simulated clusters lack. 
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Fraction of virial radius 

Figure 6. The fraction of the total bolomctric luminosity that is emitted from within the specified radius for each of the clusters 
(calculated from equation 1). Clusters with significant substructure are shown as dotted lines. 



Ponman, Cannon & Navarro (1999) argue for such preheat- 
ing by examining ROSAT observations of 25 clusters. 

The inclusion of cooling into a cosmological hydrody- 
namics simulation has proved highly successful. We have 
highlighted the necessity of cooling reasonable amounts of 
gas onto the centre of each galaxy cluster if sensible X-ray 
luminosity estimates are to be obtained. We have achieved 
this by a judicious choice of mass resolution and decoupling 
the hot and cold gas phases. A rigorous method of doing this 
in SPH is under development. We next intend to simulate 
more massive clusters in the same way, in order to extend 
our predictions into regions more accessible to observation. 
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